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ABSTRACT 

The role of feedback from massive stars is believed to be a key element in the evolu- 
tion of molecular clouds. We use high-resolution 3D SPH simulations to explore the 
dynamical effects of a single 07 star emitting ionising photons at 10^^ s~^ and located 
at the centre of a molecular cloud with mass 10^ and radius 6.4 pc; we also per- 
form comparison simulations in which the ionising star is removed. The initial internal 
structure of the cloud is characterised by its fractal dimension, which we vary between 
V = 2.0 and V = 2.8, and the roughly constant standard deviation, ct^q = 0.38, of 
its log-normal density PDF. (i) As regards star formation, in the short term ionis- 
ing feedback is positive, in the sense that star formation occurs much more quickly 
(than in the control simulations), in gas that is compressed by the high pressure of 
the ionised gas. However, in the long term ionising feedback is negative, in the sense 
that most of the cloud is dispersed with an outflow rate of up to ~ lO~^M0yr~^, on a 
timescale comparable with the sound-crossing time for the ionised gas (~ 1 to 2Myr), 
and triggered star formation is therefore limited to a few percent of the cloud's mass. 
We will describe in greater detail the statistics of the triggered star formation in a 
companion paper, (ii) As regards the morphology of the ionisation fronts (IFs) bound- 
ing the Hll region and the systematics of outflowing gas, we distinguish two regimes. 
For low V < 2.2, the initial cloud is dominated by large-scale structures, so the neu- 
tral gas tends to be swept up into a few extended coherent shells, and the ionised 
gas blows out through a few large holes between these shells; we term these Hll re- 
gions shell- dominated. Conversely, for high V > 2.6, the initial cloud is dominated by 
small-scale structures, and these are quickly overrun by the advancing IF, thereby pro- 
ducing neutral pillars protruding into the Hll region, whilst the ionised gas blows out 
through a large number of small holes between the pillars; we term these Hll regions 
pillar- dominated, (iii) As regards the injection of bulk kinetic energy, by ~ 1 Myr, 
the expansion of the Hll has delivered a mass-weighted root-mean-square velocity of 
~ 6kms~^; this represents less than 0.1% of the total energy radiated by the 07 star. 

Key words: Galaxies: ISM - ISM: nebulae - Hll regions - bubbles - Hydrodynamics 
- Stars: formation 



1 INTRODUCTION 



Infrared shells and bubbles ( [Churchwell et al. j2006 Church- 1 
1^^2008) are ubiquitous in the Galaxy and are often as- 
sociated with Hii regions, i.e. localised regions in which the 
Lyman continuum radiation (E^ > 13.6 eV) from one or 
more young, massive stars ionises the gas and heats it to 
Tj ~ 10^ K. Because the neutral gas outside the Hii region 
is much colder (T^ ^ 10 to 30 K), the Hii region is overpres- 
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sured, and expands, sweeping up the neutral gas. Molecu- 
lar clouds typically have a complicated clumpy structure, 
and so the ionising photons can penetrate to different dis- 
tances in different directions, causing highly irregular ioni- 
sation fronts (e.g. Hegmann & Kegel"2003V Consequently, 
large-scale Hii regions have complex shapes and diverse mor- 
phologies. Sometimes they appear as bright-rimmed clouds 
or perfectly round shells (Dehar- 



(Urquhart et al. 



veng et al. 2009 



2009) 



2010), other times they are dominated 



by pillar-like structures (Sugitani et al. 2002 Smith et al 



2010 Preibisch et al. 2011), or a mixture of all three fea- 
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tures ( |Koenig et al.|2008| ). These diverse morphologies must 
encode important clues as to how a massive star interacts 
with the surrounding interstellar medium (ISM). 

Feedback from young, massive stars plays a critical role 
in the evolution of molecular clouds, and of the galaxies 
in which they reside. The input of radiative and mechani- 
cal energy due to ionising radiation and stellar winds from 
massive stars is most influential during the early evolution 
of a molecular cloud, before the first supernova explodes 
after > 3Myr. In particular, the expansion of Hii regions 
may have a positive effect and trigger new star formation 



(|Elmegreen|1998||Elmegreen et al.|20 02||Whitworth Fran^ 
cis||2002 Elmegreen||2011 ), but it may also have a negative 



worth||1979] |Tenorio-Tagle|[T979l [Bodenheimer et al.||1979f 



Wilhams fc McKee|1997| [Franco et al.|199Qp Matzner (2002f , 

thereby terminating star formation, and possibly also un- 
binding the newly-formed cluster of stars ( Dale et al.|2005 ). 
Ionising feedback may also play a key role in regulating the 
star formation efficiency on galactic scales ( [McKee 1989 
Vazquez- Semadeni et al.""2010); and Hii regions and bub- 



bles excavated by ionisation may be the first step towards 
launching large scale winds and outflows from galactic disks 
( Tenorio-Tagle et al.|2003| [Hopkins et al.||2Q12 ). 

Many numerical studies concerned with positive feed- 
back, i.e. the triggering of star formation, have been per- 
formed in recent years: on individual molecular cloud cores 



(Bisbas et al. 2011 


Haworth & Harries 2012); in the con- 


text of bright 


rims and pillars (Elmegreen et al.|1995||Miao 


et al.[ [20061 


Gritschneder et al. [2009| 2010| [Ercolano & 


Gritschneder 


2011 


Mackey & Lim 2011), on parsec-scale 


molecular clouds (Mellema et al.|2006 Krumholz et al. 2007 


Arthur et al. 


2011 


Walch et al. 2011); and on whole star 



difficult to distinguish triggered from spontaneous (untrig- 
gered) star formation, even though many authors have tried 
to identify age gradients in the young stars surrounding Hii 
regions (Preibisch & Zinnecker[[2007), and to establish evi- 



dence for the Collect & Collapse scenario in swept-up shells 
(Zavagno et al.||2006 , Deharveng et al.[[2008| ^Zavagno et al 



2010). 

The negative effect of ionising feedback, i.e. its role 
is dispersing molecular clouds, has been evaluated analyt- 

has consid- 



ically by [Whitworthj ([1979f. [Matzner| f2002) 
ered cloud destruction by all stellar feedback mechanisms 
and concludes that ionising feedback will destroy normal 
giant molecular clouds in less than a crossing time. The ef- 
fect of ionsation feedback has been explored numericUy by 



Dale et al. (2005j). Based on simulations by Bonnell & Bate 
(2002), they model the impact of ionising feedback in tur- 



bulent, star-cluster- forming clouds, and find that photoioni- 
sation is able to unbind star clusters by efficiently removing 
gas quite early on during their evolution. They note that 
the mean gas density of the affected cloud is the critical pa- 
rameter. Krumholz et al. ( 2006 ) study the disruptive power 



of Hii regions using semi- analytical models, and conclude 
that lower mass clouds (M < 3x lO^Mo) may be destroyed 
in ^ lOMyr. In more massive clouds disruption by ionisa- 
tion plays a less important role ( Krumholz Matzner|2009 



Murray et al.|2010|[Fall et al.|2010t . These conclusions have 
recently been conffimed in 3D SPH simulations by |Dale| 
et al. (2012^), who study star-cluster forming clouds with 



10^ to 10^ Mo; their lower mass clouds are comparable with 
the ones we investigate in the current study. They find that 
massive clouds resist dispersal by ionising feedback, because 
it is quenched by accretion flows. Quenching of feedback by 
strong accretion is also important in the early stages of Hll 
region formation, i.e. in the ultra-compact stage (Mac Low 
et al.|2007[ [Peters et al.|2010[[201ll ). Moreover, the location 
of a massive star is important for the dispersal efficiency. 



This has first been noted by Mazurek ( 1980 ) and Yorke et al. 
(1989), and recently conffimed by Gendelev & Krumholz 
(|2012), who perform three-dimensional radiation-magneto- 
hydrodynamic simulations of blister-type Hii regions. They 
note that the recoil provided by escaping gas may also play 
a role in destroying the parental molecular cloud. 

Feedback also drives turbulent motions in molecular 
clouds. Matzner (2002) argues that Hii regions are the main 
drivers of molecular cloud turbulence, being energetically 
dominant over the combined effects of stellar winds and 
supernovae. Mellema et al. (2006) record that the level of 
molecular cloud turbulence is maintained at ^ 8kms~^ in 
their simulations of ionising feedback, but they also note 
that radial velocities account for a significant fraction of the 
total kinetic energy. [Gritschneder et al.[ ([2009 ) discuss tur- 
bulence driving in their simulations of ionising feedback, and 
report that the energetic ratio of solenoidal to compressible 
modes (Federrath et al. 2008) is significantly reduced, rela- 
tive to the fiducial ratio of 2:1. 

In this paper we simulate the development of large-scale 
Hii regions in molecular clouds having different initial fractal 



dimension V. In contrast with [Dale et al. (2005) our molec- 



ular clouds all have the same mass, radius and density PDF. 
We demonstrate that the full range of morphological features 
observed in Hii regions might be attributable to the fractal 
dimension, P, of the parental molecular cloud. As V is in- 
creased from small values P ~ 2, there is a shift from Hii re- 
gions bounded by relatively smooth extended shell-like IFs, 
with ionised gas blowing out through large holes between 
these shells, to Hii regions bounded by much more irregular 
IFs with numerous pillars protruding into the Hii region, 
and ionised gas blowing out through many small holes. The 
kinetic energy injected into the ISM is essentially indepen- 
dent of P, and the time taken to disperse the cloud is short 
(- IMyr for V = 2.8 to - 2 Myr for V = 2.0). 

The plan of the paper is the following. In section |2] we 
describe the generation of fractal molecular clouds and the 
numerical method, including the treatment of the ionising 
radiation. In section |3] we analyse the large-scale structure 
of the Hii regions obtained with different fractal dimensions, 
and in section |4] we analyse the structure of the correspond- 
ing IFs and the associated outflow rates. In section |5] we 
evaluate the injection of kinetic energy into the ISM. In sec- 
tion [6] we summarise our main conclusions. 



2 INITIAL CONDITIONS AND NUMERICAL 
METHOD 

2.1 Generation of fractal molecular clouds with 
an FFT algorithm 

It is well known that molecular clouds are rich in inter- 
nal structure, which may derive from pure turbulence (e.g. 
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Klessen|[200H |Padoan k NQrdlund|[2002t , or from gravita- 



tional amplification of noise in the density field (Bonnell 
I&: B ate 2002) or from initial turbulence leading to den- 
sity fluctuations which are subsequently enhanced by gravity 



( Heitsch et al.|2009||Ballesteros-Paredes et al.|2011|?||Walctr 



et al.||2010[|2012||Girichidis et al.||20lH [201 2^ It is unclear 
which of these processes is dominant in shaping molecular 
clouds. However, what we do know is that on scales from 
the typical prestellar core size of ^ 0.1 pc, up to ~ 100 pc 
( Bergin Tafalla|2007 ) , molecular clouds subscribe to a sta- 
tistically self-similar fractal structure (Stutzki et al. 1998 



Elmegreen Falgarone||1996| Falgarone et al.||1991| 

Apart from their applicability, we adopt fractal initial 
conditions for their practicality. Fractals with mass-size re- 



lations similar to those observed ( Larson||l98l|) are easily 
created in Fourier space ("Elme green 



2002 Shadmehri 



Elmegreen|2011 ), whilst at the same time retaining the free- 



dom to choose the fractal dimension, P, arbitrarily. In con- 
trast, when starting with a periodic box of driven or decay- 
ing turbulence (see e.g. Arthur et al.|[20TT for particularly 
nice simulations of this type), P is a consequence of the 
physics implemented, and cannot easily be adjusted. 

( 1998 ) have shown that a three- 



Stutzki et al. 



dimensional density field in which the power spectrum of 
the density fluctuations has index —n (i.e. P{k) oc has 
box-coverage fractal dimension 



along each Cartesian axis (x, y, z), where /c = 1 corresponds 
to the linear size of the box (15 pc). After generating the 
spectrum of density fluctuations in Fourier space, and trans- 
forming them to obtain ^^^^(x, y, z) in real space, this field 
is scaled according to 



p{x,y,z) = exp < 



,{x,y,z) 



(2) 



The resulting density field, p{x,y, z), has a log-normal den- 
sity PDF and a clump mass distribution in agreement with 
observations (as described in Shadmehri & Elmegreen 2011}. 
For a given spectrum of density fluctuations, changing po 
changes the standard deviation, a, of the log- normal density 
PDF, whilst leaving the underlying topology of the density 
field unchanged: small p^ produces a very wide density PDF 
(large a), and conversely large p^ produces a very narrow 
density PDF (small a). Before creating the initial SPH par- 
ticle distribution, we shift the point of maximum density to 
the centre of the computational domain and cut out a sphere 
of diameter 12.8 pc centred on this point. Finally, the box is 
partitioned with a 128^ grid, which is then populated with 
SPH particles, by placing the appropriate number of par- 
ticles randomly within each grid cell. We have checked the 
fidelity of this procedure by performing a back transforma- 
tion and recovering the correct power spectrum. 



P = 4 



(1) 



(see also [Federra th et al.||2009 ). Thus, defining V is equiva- 
lent to defining the power spectral index n, and vice versa. 
Low V corresponds to high n, i.e. a density field domi- 
nated by a small number of extended density structures; 
conversely, high V corresponds to low n, i.e. a density field 
dominated by a large number of compact density structures. 

CO line surveys are well suited to determining the struc- 
tures of molecular clouds like the ones we simulate here, be- 
cause this transition is a dominant coolant at number den- 
sities - 10^ - 10^ cm . Values of the fractal dimension es- 
timated in this way are typically P 2.3 =b 0.3 (Elmegreen^ 
&: Falgarone||1996 ). For example, observations of molecular 
clouds in the Milky Way typically yield V ^ 2.4 ([Falgarone 



et al.||1991| [Vogelaar Wakker| |1994l [Stutzki et al.||1998| 



Lee 2004; Sanchez et al. 2005), corresponding to n = 3.2. 



For Ophiuchus, Perseus, and Orion, Sanchez et al. (2007) 
obt ain a slightly higher value , P = 2.6=b0.1 (n = 2.8 ± 0.2), 
and jMivihe-Deschenes et al] ( |2010t derive V ^ 2.65 =b 0.05 
(n ^ 2.7 =b 0.1) from Herschel- SPIRE observations of the 



Polaris Flare. Sanchez et al. (2007) and Schneider et al.j 
(2011) report differences between low-mass star forming re- 



gions and massive, giant molecular clouds, based on the 
A— variance method applied to Ay maps (Ossenkopf et al.' 
2008), where Ay is the visual extinction. In M33, Sanchez 



et al. (2010) find V 2.2 to 2.5 on scales < 500 pc, and 



higher values of V on larger spatial scales. 

Here, the initial three-dimensional fractal density field 
is constructed in a periodic box, using an FFT-based algo- 
rithm ( ^Shadmehri Elmegreen 2011). The algorithm has 
three main input parameters: (i) the power spectral index, 
n = 2(4 — P), (ii) the random seed 71 used to generate differ- 
ent realisations; and (iii) the density scaling constant po (see 
Eqn. [2]below). We populate the integer modes /c = 1, .., 128 



2.2 Numerical method 



We use the SPH code Seren (Hubber et al. 2011), which 



is well-tested and has already been applied to many prob- 
lems in star formation (e.g. Walch et al.||201l" Bisbas et al. 



2011[ jStamatellos et al. 2011). We employ the standard 
SPH implementation (Monaghanj [T992t . The SPH equa- 



tions of motion are solved with a second-order leapfrog in- 
tegrator, in conjunction with an hierarchical, block time- 
stepping scheme. Gravitational forces are calculated using 
an octal spatial decomposition tree ^Barnes & Hut 1986), 
with monopole and quadrupole terms and a Gadget-style 



opening-angle criterion (Springel et al. 2001). We use the 



standard artificial viscosity prescription (Monaghan & Gin- 
gold|1983 ), moderated with a Balsara switch fBalsa ra|1995 ). 
Ionising radiation is treated with an HEALPix-based adap- 
tive ray-splitting algorithm, which allows for optimal reso- 
lution of the ionisation front in high resolution simulations 
Bisbas et ar[[2009| ). We assume the On-The-Spot ap- 



proximation and do not treat the diffusive radiation field. 

The temperature of ionised gas particles is set to 
Tj = 10^ K. The temperature of neutral gas is given by a 
barotropic equation of state. 



T{p) 



1 + 



(7-1)^ 



(3) 



where = 30 K, p^^^^ = 10~^^gcm"^, and 7 = 5/3. The 
choice of T-^ = 30 K influences the gravitational stability of 
the neutral gas swept up by the expanding Hii region. If the 
shell were able to cool to lower temperature, it would frag- 
ment more readily than currently seen in our simulations. 
For this reason we will explore a more realistic cooling func- 
tion in a future paper. With the above prescriptions for the 
temperature, and assuming solar elemental composition, the 
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ID 






^10 


^1, ^15 






g cm~3 




Myr 


P2.0/O7 


2.0 


1.29 


0.39 


0.51,0.77 


TO.2/07 


2.2 


1.17 


0.39 


0.48,0.65 


V2A/07 


2.4 


1.02 


0.38 


0.44,0.59 




2.6 


0.98 


0.38 


0.49,0.77 


V2.^/07 


2.8 


0.93 


0.37 


0.62,0.89 


TO.O/NI 


2.0 


1.29 


0.39 


0.98,0.99 




2.2 


1.17 


0.39 


0.94,0.99 


TO.4/NI 


2.4 


1.02 


0.38 


>1.0, >1.0 




2.6 


0.98 


0.38 


>1.0, >1.0 


P2.8/NI 


2.8 


0.93 


0.37 


>1.0, >1.0 



Table 1. Initial conditions. Column 1 gives the run ID (where 
the first element of the ID gives the fractal dimension, and the 
second records whether there is an ionising star, 07, or no ion- 
ising star, NI), column 2 the fractal dimension, column 3 the 
mass- weighted mean density, pM-, and column 4 the logarithmic 
standard deviation, ct^q, of the mass-weighted density PDF (see 
Fig. Column 5 gives ti and tis, where ti is the time at which 
the first sink forms, and ^15 is the time at which the fifteenth sink 
forms. 



1 L I I I I 




-23 -22 -21 -20 -19 

logio(/5/ [g cm"^]) 

Figure 1. Initial mass- weighted density PDFs for T> = 
2.0, 2.2, 2.4, 2.6 and 2.8. Ngp^^ is the number of SPH particles in 
each density bin and N^q^^ is the total particle number. Means 
and standard deviations are given in Table [l] 



isothermal sound speed is Cj = 12kms~^ in ionised gas, and 
— 0.32 km s~"^ in neutral gas at low density {p <^ Pcrit)- 
Although we are not concerned with triggered star for- 
mation in this paper, we introduce sinks at density peaks 
above yOgi^K — 10~^^ gcm~^, provided that the density peak 
in question is at the bottom of its local gravitational poten- 



tial well (see Hubber et al. 2012[ in preparation). Since 



PsiNK ^ PcRiT ' ^ condensation that is converted into a sink 
is always well into its Kelvin-Helmholtz contraction phase. 
Once formed, a sink is able to accrete gas smoothly from its 
surroundings and thereby grow in mass. 



2.3 Initial conditions 

We consider a single cloud mass, 10^ M©, and a single cloud 
radius, 6.4 pc, so the volume- weighted mean density of the 
cloud is pv = 0.62 x 10~^^ gcm""^, the escape velocity from 



the surface of the cloud is 
freefall time is tp 
a total number of N^^ 



:3.7kms , and the mean 
3 Myr. Every cloud is modelled with 
2.5 X 10^ SPH particles, re- 



sulting in a mass resolution of O.4M0 (Bate Burkert 



1997). The only parameter we vary is the initial fractal di- 
mension, V = 2.0, 2.2, 2.4, 2.6 and 2.8. For each value of 
P, we generate a single realisation, using the same random 
seed, 7^= 1, so that the density fields in the five different 
clouds have the same underlying topology. We also adjust po 
so that the (approximately log-normal) density PDFs of the 
five clouds have roughly the same standard deviation (see 
Fig[l]for the initial mass- weighted density PDFs). We then 
simulate the evolution of each cloud twice, once with a single 
07 star at the centre of the cloud, emitting ionising photons 
at A/'Lyc — 10^^ s~^, and once with no ionising star. Table 
[1] lists the ID of each simulation, the fractal dimension, the 
mass- weighted mean density, pM, the standard deviation of 
the logarithmic density PDF, ct^q, and the times, and t^^, 
at which the first and fifteenth sinks are created (a measure 
of how rapidly star formation occurs) . We note that pM is in- 
evitably larger (up to a factor of 2) than pv - In the lefthand 



column of Fig. |2] we show initial column-density images of 
all five clouds, as seen projected on the (x,y)-plane. We note 
that as V is increased, the scale of the dominant features 
in these images becomes smaller. In addition to these ten 
simulations involving lumpy clouds, we perform a control 
simulation with a uniform-density cloud of the same mass 
and radius: this has the ID 'control'. 



3 GENERAL EVOLUTION 
3.1 Morphology 

The morphology of the evolving Hii region is strongly de- 
pendent on the fractal dimension of the initial molecular 
cloud. In the two central columns of Fig. [2] we show column- 
density images of the five clouds as projected on the (x,y)- 
and (x,z)-planes, at time t — 0.66 Myr. By this time there 
has been considerable interaction between the Hil region and 
the molecular gas, and in all cases star formation has been 
triggered (turquoise dots mark sink particles); we will anal- 
yse the triggering of star formation in a companion paper 
(Walch et al., in prep.). 

Although the simulated Hii regions have highly chaotic 
structures (which could be explored further by invoking dif- 
ferent random seeds, 7^), there are also discernable statisti- 
cal trends with changing V. For low T>, the dominant density 
structures in the initial cloud have large coherence lengths, 
and consequently the Hii regions tend to be bounded by 
extended, relatively smooth shells, and their external ap- 
pearance can be strongly dependent on viewing angle. Be- 
tween the shells are large holes, through which the ionised 
gas can escape into the ISM. We term these Hii regions 
shell- dominated. Conversely, for large P, the dominant den- 
sity structures in the initial cloud have relatively short co- 
herence lengths, and consequently the Hll regions tend to 
have a much more irregular boundary, with many neutral 
pillars protruding into the ionised gas. There are many 
small holes through which the ionised gas escapes. We term 
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these Hii regions pillar- dominated. The transition between 
shell-dominated and pillar-dominated Hii regions occurs at 
V ~ 2.4, which is the average fractal dimension of many 
observed molecular clouds 



2.1) 



(see section 

In the righthand column of Fig. [2] we show, for com- 
parison, column-density images projected on the (x,y)-plane 
for the simulations performed without an ionising star. Un- 
der this circumstance the evolution proceeds on a free-fall 
timescale, ^ 3Myr, so that after 0.66 Myr there has been 
much less change to the density field. 

We conclude that the main requirement for the forma- 
tion of pillars is the existence of small-scale structure in the 
density field into which the IF advances. As V increases, 
and therefore n decreases, there is more small-scale struc- 
ture and more pillars are formed. Since, for turbulent gas 
with a fixed ratio of solenoidal to compressive modes, small- 
scale structure should become more abundant with increas- 
ing Mach Number, this conclusion agrees with the findings 



of Gritschneder et al. (2010). They have performed simula- 



tions of the ionisation of a slab of turbulent interstellar gas, 
and invoke different Mach numbers (where the Mach Num- 
ber is measured relative to the sound speed in the neutral 
gas). They conclude (i) that the most favourable regime for 
the formation of pillar-like structures is Mach 4 to 10, and 
(ii) that the typical size of a pillar is set by the scale of the 
largest turbulent perturbation. 



For the fractal clouds, R^^ is initially smaller than , 
because the 07 star is placed at the densest point in the 
initial density field. However, as the IF advances away from 
this point, the density falls below pv, so that R^^ expands 
faster than oc t^^^ and catches up with the uniform-density 
solution after ^ 0.2 Myr. Thereafter R^^ tracks the uniform- 
density solution quite closely (±0.5 pc) for all V. For fractal 
clouds, there are two competing effects: (i) R^^ is reduced 
because the volume- weighted mean square density, (p^)^, 
is higher than for the uniform density in the CONTROL sim- 
ulation, and therefore the ionising photons should get used 
up balancing recombination in a smaller volume; (ii) is 
increased because the IF can break out of the cloud more 
readily, through paths of low emission measure. Effect (i) ap- 
pears to dominate for low T>, making R^^ smaller; and effect 
(ii) appears to dominate for high P, making R^^ larger. 



4 STRUCTURAL ANALYSIS OF SHELLS 

From Fig. |2] it is evident that, over and above stochastic 
variations, there is a systematic progression in the struc- 
tural properties of the Hii regions with increasing V. In this 
section, we analyse the shell structures more quantitatively, 
and compute outflow rates of hot and cold gas, in order to 
derive cloud destruction timescales. 



3.2 Radius of the Hii region 

The mean radius of the ionisation front, R^^ (see Fig. [3| 



an observational quantity which is easily obtained ( Church- 
well et ar]|2006| , and frequently used in statistical studies 
(e.g. in stud ies of triggered star formation by Thompson 
et al. 2011). Fig. [s] shows the evolution of R^^ obtained 
with different T>. i^jp is estimated by averaging the radii 
of all SPH particles having temperature between 500 K and 
5000 K (see |Bisbas et al.|2009| . 

Given a uniform-density ambient medium, as in the 
CONTROL simulation, the IF advances to the Str0mgren ra- 
dius. 



3A/;.c < 

4:7V a^X'^ p2 



1/3 



in a few recombination times, 



ttg X pv 



0.4 kyr. 



(4) 



(5) 



Here = 1.7 x 10 

2.6 X 10~^^cm^s"^ 



g is the mass of an hydrogen atom, 
is the Case B recombination coef- 
ficient (i.e. invoking the On-The-Spot Approximation, and 
therefore neglecting recombinations straight into the ground 
state), and X = 0.7 is the fraction by mass of hydrogen. For 
the CONTROL simulation we derive R^ — 2.3 pc. Thereafter 
the IF quickly switches to D-type, and the radius is given 
approximately by ( |Shu|1991 ) 

1 + i> ^ • (6) 



4Rj 



This solution is very well reproduced in the control simu- 
lation, until the IF reaches the edge of the cloud (see dashed 
and dotted lines on Fig.js]). 



4.1 Shell structure 

Fig. [4] shows Hammer projections of the column-density, as 
seen from the ionising star, at t = 0.66 Myr. If we adopt 
spherical polar coordinates, (r, 0, 0), centred on the ionising 
star, then the false colour on the Hammer projection repre- 
sents the column-density, S, as seen in each direction, {0, 0), 
from the ionising star. 



E(0,</.) 



p(r, 0, (j)) dr . 



(7) 



For low P, the Hammer projection is dominated by a small 
number of extended column-density features, which corre- 
spond to the extended convex shell-like segments of the ion- 
isation front seen in the top rows of Fig. [2] because there 
is only a small number of extended features these Hii re- 
gions can look quite different to different external observers 
(compare the middle images on the top two rows of|2|. Con- 
versely, for high P, the Hammer projection is dominated by 
a larger number of more compact column-density features, 
corresponding to small concave bright-rimmed clouds and 
pillars, as seen in the bottom rows of Fig. [2] because there 
is a large number of compact features, these Hii regions look 
rather similar to different external observers (compare the 
middle images on the bottom two rows of Fig. [2|. 

The white contours on Fig. |2] demark the separatrices 
between directions in which the Hii region is ionisation- 
bounded (i.e. rays along which the ionising photons run 
out before reaching the edge of the cloud), and directions 
in which the Hii region is density-bounded (i.e. rays along 
which there is still ionising radiation left at the edge of the 
computational domain). Table[2]gives the fraction of the sky. 
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Figure 2. From top to bottom, the different rows correspond to T> = 2.0, 2.2, 2.4, 2.6 and 2.8. The lefthand column shows column- 
density images of the initial fractal clouds, projected onto the (x, 2/)-plane. The middle two columns show column-density images at 
t = 0.66 Myr from the simulations with ionising radiation, projected onto the (x,y)- and (x, 2;)-planes. The righthand column shows 
column-density images at t = 0.66 Myr from the simulations without ionising radiation. The colour-table on the right gives the logarithm 
of column-density in Mq pc~^. 
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0.4 0.6 
Time [Myr] 

Figure 3. The solid curves give the mean radius of the ionisation 
front, Rip, as a function of time, for different V (see key for 
colour coding of different V). The dashed curve represents the 
uniform-density CONTROL simulation, and the dotted curve shows 
R^p oc f^/^. The thin, dashed horizontal line marks the Str0mgren 
radius of the CONTROL run, and the thin, dotted horizontal line 
marks the initial cloud radius. 





V = 


= 2.0 


V = 


= 2.8 


t = 0.30 Myr 


/iB " 


= 0.72 


flB " 


= 0.80 


t = 0.66 Myr 


/iB " 


= 0.62 




= 0.45 



Table 2. The fraction of the sky that is ionisation bounded, as 
seen from the ionising star, /jg, for V = 2.0 and 2.8, at times 
t = 0.30 and 0.66 Myr. 



seen from the ionising star, that is ionisation bounded, 
fu 



= 1- [ 

47T / 

TPt 



(8) 



for V = 2.0 and 2.8, at times t = 0.30 and 0.66 Myr. The 
large-scale holes that dominate the low V cases form rel- 
atively quickly after the ionising star switches on, as their 
formation is based on the pre-existence of coherent regions 
of low density which are too extended to be filled in sig- 
nificantly by the erosive flows at their boundaries. However, 
because the net length of these boundaries is relatively short, 
the net rate at which the remaining neutral gas is eroded is 
slow. Consequently /jg starts small but decreases slowly in 
these cases. In contrast, the small-scale holes that dominate 
the high V cases are initially choked because they are small 
and therefore the erosive flows at their boundaries are much 
more effective at inhibiting the escape of ionising radiation. 
However, because the net length of these boundaries is rel- 
atively large, the net rate of erosion is high. Consequently 
/jg starts high, but decreases rapidly in these cases. 



4.2 Outflow rates 

Fig. [6] shows the total rate at which mass flows out through 
the initial boundary of the cloud (at radius 6.4 pc), as well 
as the rate for the ionised gas only, both as functions of 



lO^M0Myr"\ and 10 to 25% of this is ionised gas. The 
dependence of M on P is rather weak. 

Fig. [7|l shows the root-mean-square velocity for all the 
gas and for the neutral gas only, plus the mean radial ve- 
locity of the neutral gas, all as functions of time, t, for the 
runs with an ionising star at the centre. By ^ 1 Myr, the 
radial velocity of the neutral gas is approaching the escape 
speed from the cloud surface, and the ionised gas is travel- 
ling outwards even faster. We conclude that the cloud will 
be dispersed by ^ 1 to 2 Myr, (where the longer time corre- 
sponds to lower V). Fig. [7^ also shows the root-mean-square 
velocity for all the gas when there is no ionising source. In 
this case the velocities are generated by self-gravitational 
acceleration, and are much smaller. Therefore we can affirm 
that it is the action of the ionising star that dominates the 
cloud dynamics. 



5 TURBULENCE DRIVING & DENSITY 
EVOLUTION 

5.1 Driving turbulence 

The injection of bulk kinetic energy by expanding Hii regions 
may make a significant contribution to the driving of tur- 
bulence in molecular clouds. Gritschneder et al.| ( |20Q9| ) have 
shown that ionising radiation is effective in driving turbu- 
lence in the surrounding ISM; in particular they find that 
the fraction of turbulent energy contained in compressional 
modes is enhanced relative to the equipartition value of 1/3. 
The simulations presented here admit an even more unequiv- 
ocal opportunity to evaluate the role of ionising radiation in 
driving turbulence, since there are no initial turbulent ve- 
locities, and therefore we do not need to take account of 
dissipation. The total kinetic energy must be attributed to 
ionising feedback or to self-gravity. By performing compar- 
ison simulations with no ionising source, we can estimate 
accurately the contribution from ionising feedback. 

From Fig. |7|i we see that, to a good approximation, 
the root mean square velocity increases linearly with time, 
both with (full curves) and without (dashed curves) ionis- 
ing feedback. However, the rate without ionising feedback, 
^ lkms~^ Myr~^, is much smaller than the rate with ion- 
ising feedback, ^ 6kms~^Myr~^. The corresponding rates 
of injection of kinetic energy, expressed per hydrogen nu- 
cleus, are 16,^tm 4 x 10~^^ ergs~^ H~^ without ionising 



feedback, and ii^ 



1.4 X 10 ^^ergs ^H ^ with ionising 



feedback; for comparison, the heating rate due to ionisation 
is - 5 X 10"^^ergs-^H-^ (p/6.2 x 10"^^gcm-^) - 

where we have normalised the density to the mean value in 
the initial cloud. 

Fig.|7|i also shows the root mean square velocity (dotted 
curves) and the radial velocity (dash-dot curves) for the neu- 
tral gas only, in the case where there is ionising feedback. We 
see that even the neutral gas acquires a large amount of bulk 
kinetic energy; and that a large part of this is invested in or- 
dered expansion. The radial velocity is approaching 4 km s~^ 
at the end of the simulation, which is comparable to the 
bubble expansion velocity found in IC 1396, a shell-like HII 



time. By 1 Myr, the total rate is 



0.5 to 1.5 X 



region with a comparable radius (Patel et al. 1995). How- 
ever, even if we discount the part invested in ordered expan- 
sion, the remaining random kinetic energy ("turbulence") is 
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Figure 4. Hammer projections of the column-density, ^{0,(j)), 
as seen from the ionising star at time t = 0.66 Myr. From top 
to bottom, V = 2.0, 2.2, 2.4, 2.6 and 2.8. The white contours 
mark the separatrices between ionisation-bounded directions and 
density-bounded directions. 
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Figure 5. The fraction of the sky, as seen from the ionising star, 
that is ionisation-bounded, /jg , as a function of time, t. From top 
to bottom, V = 2.0, 2.2, 2.4, 2.6 and 2.8. 




0.0 0.2 0.4 0.6 0.8 1.0 

Time [Myr] 



Figure 6. Mass outflow rates through the initial surface of the 
sphere at radius 6.4 pc for V = 2.0, 2.2, 2.4, 2.6 and 2.8. Solid 
lines show the total outflow rate, whereas dashed lines show the 
amount of outflowing ionised gas. 



injected into the neutral gas at a greater rate by ionising 
feedback (^-^ 2 — 4 km s""*^), than by gravity when there is no 
ionising feedback. Moreover, after ^ IMyr the mean radial 
velocity of the neutral gas is approaching the escape veloc- 
ity from the initial surface of the cloud (3.7km s~^), and 
therefore the cloud is being effectively dispersed. 

We note that some of the kinetic energy invested in or- 
dered expansion could be transformed into turbulent energy 
when the expanding shell interacts with the surrounding in- 
terstellar medium. For instance, in their three-dimensional 
radiation- hydrodynamics simulations, studying the ionisa- 
tion of turbulent clouds without self-gravity, [Mellema et al.| 
(2006| derive root-mean square velocities for the cold gas. 
In particular, they find that, after the initial ^ 10^ yr, the 
difference between the one-dimensional rms velocity of the 
cold gas and the average expansion velocity of the cold gas is 
roughly constant at ^ 4kms~^. Thus, the pre-existing tur- 
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bulence of 2 km seems to be maintained and enhanced 
by the ionising feedback. Therefore, the amount of energy, 
which is found to be transformed into random motions is 
comparable to our simulations. 

Fig. [T}) shows the evolution of the bulk kinetic energy, 
both with (full curves) and without (dashed curves) ionis- 
ing feedback; also plotted is the net energy input in the form 
of ionising radiation (dotted curve, scaled down by 10~^). 
By ~ 1 Myr, the total radiative energy injected by the ion- 
ising source is 7 x 10^^ erg. A part of the injected energy 
is transformed into kinetic energy ^ 3 x 10^^ erg. In the 
run without ionisation, only gravitational potential energy 
may be transformed into kinetic energy. We find that this 
results in ^ 10^^ erg. Thus ionising feedback is much more 
effective than gravity at generating bulk kinetic energy, at 
least on this timescale (^ IMyr), but still very inefficient 
0.05%). Similar results have been found by Dale et al. 
( |2005 ), who model somewhat smaller, less massive clouds, 



having comparable mean densities; they do not place an ion- 
ising source at the center, but instead evolve a molecular 
cloud with turbulent initial conditions to form a massive star 
self-consistently. In radiation-magneto-hydrodynamic simu- 



lations Gendelev & Krumholz (2012) have recently shown 
that the inclusion of a magnetic field may enhance the in- 
put of kinetic energy to the surrounding clouds, especially if 
the ionising star forms at the edge of a molecular cloud and 
generates a blister- type Hii region. However, the energy only 
increases by a factor of ^ 2 in this case, so the efficiency is 
still < 0.1%. 



5.2 Timescales for spontaneous star formation 

Since no star formation occurs in the comparison simulations 
during the first ^ 0.9 Myr, we can infer that, up to this point, 
all star formation in the simulations with ionising feedback 
is triggered. To understand this result, consider the free-fall 
time in gas with density p. 



32Gp 
2 Myr 



1/2 



10 ^^gcm 



-1/2 



(9) 



and the sound crossing time for a cloud of neutral gas with 
radius R, 



R 



3 Myr 



(10) 
(11) 



Spontaneous (i.e. untriggered) star formation cannot oc- 
cur before 0.9 Myr, unless there are coherent lumps of gas 

The first condition 
and 

the second condition requires coherent lumps with radius 
R > 0.3pc(p/5 X 10~^^gcm"^)"^/^. From inspection of 
the density PDFs in Fig. [l] it is clear that there is very lit- 
tle gas at such high densities, and what there is is unlikely 
to be in sufficiently large coherent lumps. 



with tpp < 0.9 Myr and tpp ^ ^sc,n 
requires lumps with density p > 5 x 10~^^gcm~^ 



5.3 Density evolution 

Density PDFs can be used to characterise the structure 
of molecular clouds (e.g. [Padoan k Nordlund||2002| ). Figjs] 
demonstrates the influence of ionising feedback on the den- 
sity PDFs of the molecular clouds with T> — 2.0 and V — 2.8. 

Without ionising feedback the evolution is very moder- 
ate, and simply involves a gradual broadening of the PDF, 
until the time approaches ^0.9 Myr. At this juncture a high- 
density tail suddenly develops, for V = 2.0, signalling the 
onset of collapse and spontaneous star formation. This is 



et al. 2009 2011 Lombardi et al. 2010) and theoretical (Krit- 



in good agreement with recent observational ( [Kainulainen 



suk et al.||2011 IGirichidis et al.|2012) findings. 



With ionising feedback, the evolution of the density 
PDF is much more rapid and distinctive. A narrow low- 
density peak develops at densities p < 10~^^gcm~^, repre- 
senting an expansion wave driven by hot ionised gas escaping 
into the surrounding ISM. At the same time a much broader 
high-density peak forms at densities p ^ 3 x 10~^^ gcm~^, 
corresponding to neutral gas compressed by the overpres- 
sure of the Hii region; this neutral gas is the material for 
triggered star formation. 



5.4 Mass evolution 

Fig. [9] shows the total mass in very low-density gas, M^q^, 
and in very high-density gas, M-^^^-^ , as a function of time, 
for T> = 2.0, 2.2, 2.4, 2.6 and 2.8. Here very low density 
means the escaping ionised gas in the expansion wave with 
p < 10~^^gcm~^, and very high density means gas which 
is sufficiently dense to couple thermally to the dust, p > 
6 X 10~^^gcm~^, since this gas is very likely to condense 
into new stars. 

Without ionising feedback, the mass in low-density gas 
increases with time as the density substructure relaxes, es- 
pecially for clouds with high V. After ^ 1 Myr, M^q^ ap- 
proaches 100 M© for all V. M^^^^ is negligible in all runs 
until t > 0.9 Myr, when, for low V < 2.4, it increases rapidly, 
signalling the onset of spontaneous star formation. 

With ionising feedback M^qw ^i^^^ even more rapidly 
than without ionising feedback, as gas is blown out into 
the surrounding ISM, and the cloud is dispersed. By t = 
0.66 Myr, M^^^ has already reached 500 M©, and by t > 
0.85 Myr, M^q^ has passed 10^ M©. At the same time, 
^HiGH increases rapidly up to ^ 0.4 Myr, and then satu- 
rates at ^ 100 M0, for all T>. M^^^-^ is the mass that is 
available to form stars during the next 0.1 Myr; note that 
it does not include the mass that is already in sinks. 

Thus, although ionising feedback triggers star forma- 
tion, it also disperses the cloud on a comparable timescale. 
Consequently, triggered star formation only accounts for a 
few percent of the cloud's mass. Any other star formation 
must occur spontaneously, promptly after the ionising star 
switches on. 



6 CONCLUSIONS 

We have used high-resolution 3D SPH simulations to explore 
the effect of a single 07 star emitting photons at 10^^ s~^ 
and located at the centre of a molecular cloud with mass 
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Figure 7. Left panel: The full and dashed curves show the root-mean-squared velocity for all the gas in, respectively, the simulations 
with and without ionising feedback. The dotted and dot-dash curves show, respectively, the rms velocity and the radial velocity, for the 
neutral gas only in the simulation with ionising feedback. The thin horizontal line demarks the escape velocity from the surface of the 
initial cloud. Right panel: The full and dashed curves show the total kinetic energy for all the gas in, respectively, the simulations with 
and without ionising feedback. The dash-dot line shows the total energy input in the form of ionising radiation, scaled down by 10~^. 



10^ Mq and radius 6.4 pc. We focus on the shell structure 
and the dynamical impact of the ionising radiation on the 
surrounding molecular cloud as a function of its fractal di- 
mension, T>. We find that, although the shell morphology is 
strongly dependent on P, global parameters like the total 
outflow rate, the mass in high and low density gas, the in- 
jected kinetic energy, and the average bubble radius are in 
general almost independent of T>. 

Concerning the structures of the evolving Hii regions, 
there is a morphological transition as V is increased from 
low to high values. Clouds with low P < 2.2 are dominated 
by extended density structures. Consequently, when the ion- 
ising star switches on they develop (i) large holes, through 
which ionised gas escapes into interstellar space, and (ii) be- 
tween the holes, extended, rather smooth, shell- like struc- 
tures where the ionisation front is bound by high density 
gas. We call these clouds shell- dominated. Conversely, clouds 
with high P > 2.6 are dominated by compact density struc- 
tures. Consequently, when the ionising star switches on, they 
develop (i) a rather uniform network of small holes through 
which the ionised gas can escape into the surrounding in- 
terstellar medium, and (ii) between these small holes there 
are numerous pillar-like structures protruding into the Hii 
region. Therefore, we call these clouds pillar- dominated. The 
transition between shell-dominated and pillar-dominated 
Hii regions occurs for clouds with fractal dimension V ^ 2.4, 
which is close to the mean observed fractal dimension in 
Galactic molecular clouds. The Hii region forming in a cloud 
with D = 2.4 nicely resembles observed Hii regions (Walch 



et al.|[20TT ). We note that the presence of pillars is a sign 
of there having been small scale density structures in the 
parental molecular cloud; these might, for example, have 
been generated by turbulence. 

We discuss the influence of ionising feedback on the ki- 
netic energy budget, and on the density structure of the 
clouds. Overall, the conversion of the energy of ionising ra- 
diation into kinetic energy is extremely inefficient, < 0.1%. 



However, on the short timescale before the cloud is dis- 
persed, ^ 1 to 2Myr with M^q^ ^ lO^^Moyr"^, ionising 
feedback is much more effective at injecting kinetic energy 
than self-gravity; in comparison, for simulations without ion- 
ising feedback, the amount of kinetic energy injected after 
^ 1 Myr is 40 times smaller. The kinetic energy injected 
by ionising feedback should not be termed turbulent, since 
much of it is invested in the ordered expansion that disperses 
the cloud. 

Thus, ionising feedback is very effective at triggering 
star formation, in the sense of accelerating the rate of star 
formation, but it also terminates star formation by blowing 
the cloud apart, on a timescale ^ 1 to 2 Myr. Consequently, 
triggered star formation is rather short-lived, and only con- 
verts a few percent of the cloud's mass into stars, thus limit- 
ing the total star formation efficiency of a molecular cloud. 
Any other young stars formed in the cloud must therefore 
have been formed before or promptly after the ionising feed- 
back switched on. 
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